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various  two-dimensional  problems  of  the  dynamic  loading  of  a  slab  are  splved 
for  a  material  characterization  that  is  elastic-viscoplastic  and  exhibits 
anisotropic  work-hardening.  The  governing  constitutive  equations  are  based  on 
a  unified  formulation  which  requires  neither  a  yield  criterion  nor  loading  or 
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Induced  by  the  plastic  deformation  history.  The  theory  also  considers  plastic 
compressibility  which  depends  on  the  extent  of  anisotropy.  A  numerical 
procedure  for  solving  equations  is  developed  which  Incorporates  the  history 
dependent  anisotropic  hardening  effects.  Cases  considered  are  the  dynamic 
penetration  of  a  slab  by  a  rigid  cylindrical  indenter,  and  a  distributed  force 
rapidly  applied  over  part  of  the  slab  surface.  Both  conditions  of  fixed  and  free 
rear  surfaces  of  the  slab  are  examined.  A  uniaxial  problem  is  also  considered 
in  which  different  bases  for  the  anisotropic  hardening  law  are  examined. 
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DYNAMIC  RESPONSE  OF  A  SLAB  OF  ELASTIC-VISCOPLASTIC 
MATERIAL  THAT  EXHIBITS  INDUCED  PLASTIC  ANISOTROPY 


J. 


by 


Aboudi^ 


and  S.R. 


ABSTRACT 


Various  two-dimensional  problems  of  the  dynamic  loading 
of  a  slab  are  solved  for  a  material  characterization  that  is 
elastic-viscoplastic  and  exhibits  anisotropic  work-hardening. 
The  governing  constitutive  equations  are  based  on  a  unified 
formulation  which  requires  neither  a  yield  criterion  nor 
loading  or  unloading  conditions.  They  include  multi¬ 
dimensional  anisotropic  effects  induced  by  the  plastic  de¬ 
formation  history.  The  theory  also  considers  plastic  compress¬ 
ibility  which  depends  on  the  extent  of  anisotropy.  A  numerical 
procedure  for  solving  the  equations  is  developed  which  in¬ 
corporates  the  history  dependent  anisotropic  hardening  effects. 
Cases  considered  are  the  dynamic  penetration  of  a  slab  by  a 
rigid  cylindrical  indenter,  and  a  distributed  force  rapidly 
applied  over  part  of  the  slab  surface.  Both  conditions  of 
fixed  and  free  rear  surfaces  of  the  slab  are  examined.  A 
uniaxial  problem  is  also  considered  in  which  different  bases 
for  the  anisotropic  hardening  law  are  examined. 


1  Associate  Professor,  Department  of  Solid  Mechanics,  Materials 
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INTRODUCTION 


A  unified  theory  for  elastic-viscoplastic  isotropic 
work-hardening  materials  which  has  the  property  that  it  requires 
neither  a  yield  criterion  nor  loading  or  unloading  conditions 
has  been  proposed  by  Bodner  and  Partom  [1].  In  this  formulation, 
both  elastic  and  inelastic  deformations  are  present  at  all 
stages  of  loading  and  unloading.  This  theory  was  implemented 
in  several  two-dimensional  dynamic  indentation  problems  by 
Aboudi  [2]  and  extensive  investigations  were  made  concerning 
the  effect  of  the  viscoplastic  mechanism. 

The  isotropic  theory  ignores,  however,  the  existence  of 
the  Bauschinger  effect  which  expresses  the  observation  that 
material  hardening  in  a  given  direction  of  stressing  causes  a 
reduction  of  hardening  in  the  opposite  direction.  In  a  multi- 
axial  situation,  the  material  develops  changes  in  resistance 
to  plastic  flow  (hardening)  which  vary  in  orientation  as  well 
as  with  the  sign  of  the  stresses  (or  the  inelastic  deformation 
rates) .  The  resistance  to  plastic  flow  is  therefore  anisotropic 
and  depends  upon  the  complete  loading  history. 

In  order  to  incorporate  these  effects,  the  theory  was 
generalized  by  Stouffer  and  Bodner  [3]  to  include  anisotropic 
work-hardening.  This  is  achieved  by  the  development  of  suitable 
evolutionary  equations  for  the  inelastic  state  variables  which 
represent  the  varying  resistance  of  the  material  to  plastic 
flow  in  direction  and  in  orientation.  According  to  the  equations 
of  [3] ,  part  of  the  rate  of  change  of  a  state  variable  is 


isotropic  and  the  remaining  part  depends  upon  the 
relative  sign  and  orientation  of  the  current  rate  of 
deformation  vector.  As  a  consequence  of  the  induced 
anisotropic  hardening,  plastic  flow  of  an  initially  iso¬ 
tropic  material  becomes  anisotropic  and  compressible. 

In  the  present  paper,  the  theory  of  anisotropic 
hardening  is  implemented  on  two-dimensional  dynamic  problems 
to  investigate  the  effect  of  non-isotropic  work -hardening. 

A  numerical  procedure  is  proposed  which  generalizes  the 
previous  one  described  in  [2]  by  including  the  dependence 
of  hardening  on  the  complete  history  of  deformation.  The 
method  is  applied  to  a  class  of  dynamic  indentation  problems 
in  which  the  contact  region  between  the  indenter  (assumed 
to  be  rigid)  and  the  elastic-viscoplastic  solid  is  not 
known  in  advance  but  must  be  determined  as  part  of  the 
solution.  The  proposed  procedure  is  also  applied  to  a 
thick  slab  which  is  excited  by  an  extended  time -dependent 
normal  loading  on  its  surface,  while  the  other  surface  is 
either  rigidly  clamped  or  traction-free.  The  effects  of  the 
induced  anisotropic  hardening  and  the  consequent  anisotropic 
plastic  flow  are  shown  in  each  case,  and  comparisons 
are  made  with  corresponding  results  obtained  from  the 
isotropic  theory. 

Finally,  the  method  is  illustrated  in  a  uniaxial 
problem  of  am  elastic-viscoplastic  rod  which  is  impacted  by 
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a  time -dependent  loading.  Here  we  also  investigate  the 
effect  of  different  reference  bases  for  the  non-isotropic 
hardening  law,  i.e.  the  physical  quantities  that  could 
govern  hardening.  These  could  be  either  the  current  plastic 
deformation  rate  (as  in  [3]),  or  the  current  direct  stress. 
Compressibility  effects  associated  with  each  case  are  also 
calculated,  and  the  results  are  compared  with  each  other  and 
with  those  of  the  isotropic  theory. 
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BASIC  THEORY 


The  constitutive  equations  for  elastic-viscoplastic 
materials  with  isotropic  work-hardening  proposed  by  Bodner 
and  Partom  [ ] ]  have  the  basic  property  that  no  yield  criterion 
nor  loading  or  unloading  conditions  are  required.  In  the 
framework  of  these  equations,  both  elastic  and  inelastic 
deformations  are  present  at  all  stages  of  loading  and  unloading 
although  one  type  of  deformation  might  be  negligibly  small  in 
some  loading  regines.  For  infinitisimal  strains  (con¬ 

sidered  in  the  present  case) ,  the  constitutive  equations  of  the 
material  can  be  described  by  separating  the  total  strain  rate 
components  into  elastic  (reversible)  and  plastic  (non-re vers ible) 
strain  rates  as  follows 

*ij  =  *ij)+  i,j  =  1,2,3  U) 

where  the  strains  are  given  by  e. .  =  (u.  .  +  u.  . ) /2  with  u. 

1 » D  i 

being  the  components  of  the  displacement  vector  and  the  dots  _ 
representing  time  derivatives.  «•« 

The  elastic  strain  rates  £^®^are  related  to  the  stress 
rates  according  to  the  usual  Hooke's  Law  (assuming  elastic 
isotropy) , 

-  *ij/2y  -  X6kk6ij/[2y(3X+2u)]  (2) 

where  X,  y  are  the  Lam£  constants  of  the  material  and  6^  is 
the  Kronecker  delta. 


JfcM'NMmin** 
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In  the  most  general  anisotropic  plastic  flow  rule,  the 
plastic  strain  rates  are  related  to  the  stresses  according 
to  [3], 


P)  =  a.  s 
ei j  AijkABk* 


i » j  » k ,  £  -  1,2,3 


where  is  a  fourth  order  tensor  valued  function  and 

■ij  denotes  the  deviators  of  the  stress  tensor,  i.e., 
sij  =  0ij  “  akk'Sij/^3'  The  synunetry  properties  of  and 

Sj^j  imply  that 


Aijki  Ajik&  ^ijik 


Following  Stouffer  and  Bodner  [3] ,  let  us  rewrite  (3)  in  a 
6 -dimensional  vector  space  in  the  following  form  (the  summation 
convention  does  not  apply  to  this  paper  on  Greek  letters) : 


where 


Da 


U1  ell 

°2  e22 

D  =  £(p) 

°3  e33 

d4  -  /7 

°5  "  *  «iS* 
»6  -  *  ^ 


ct  /  3-1  /  •  •  #  /  6 


s~  .j  , 


and  Aag  is  a  6  x  6  tensor. 


in  the  6-dimensional  space,  whose  bases  are  e  ,  in  which  the 


to  the  Prandtl-Reuss  equations  for  isotropic  plastic  flow,  (7) 
implies  that  the  plastic  deformation  rate  components  Da  are 
parallel  to  the  corresponding  de viator  T  .  On  the  other  hand, 


the  anisotropic  plastic  flow  rule  (7)  does  not  lead  to 


A  which  implies 


isotropic  plastic  flow,  however,  A 


A  unified  elastic-viscoplastic  non-isotropic  work¬ 
hardening  theory  is  achieved  by  assuming  that  each  non-zero 
component  A^  (i.e.  the  diagonal  terms  Aaa)  can  be  expressed 
in  terms  of  a  single  valued  scalar  function  of  the  second 


invariant  of  the  deviatoric  stress  tensor  J~ (J 


Zag  =  0  for  a  M  so  that  each  non-zero  component  ZQg 
directly  to  a  non-zero  component  A  ft,  i.e. 
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aa 


=  F(J0,T,Z  ) 
2  aa 


(8) 


The  quantities  Z  ^  are  referred  to  as  the  hardness  variables 
of  the  material  which  express  its  resistance  to  plastic  flow. 
They  are  used  to  describe  the  state  of  the  material  micro¬ 
structure  at  any  time  t  and  thereby  depend  upon  the  complete 
history  of  the  deformation  up  to  the  current  time. 


Due  to  the  Bauschinger  effect,  the  response  of  the 
material  in  tension  and  compression  or  shear  to  the  left  and 
right  are  different  as  a  result  of  which  the  material  develops 
anisotropic  plastic  flow.  Accordingly,  we  define 


Z 

aa 


(t) 


=  < 


Zaa(t) 


Zaa(t) 


Da(t)  >  0 
Da(t)  <  0 


(9) 


where  Z+  (t)  and  Z~  (t)  are  the  hardness  variables  for  positive 

CXvJC  CXCX 

and  negative  plastic  deformation  rates,  respectively,  ( following  [3]). 


A  specific  form  for  the  evolution  laws  for  ZaQ((t)  is 
given  in  [3]  in  the  form 


Z  (t) 
aa 


q2(t)  +  (l-q)z  ra(t) 


qi(t)  -  (l-q)z  ra(t) 


Da(t)  >  0 


Da(t)  <  0 


(10) 


where  z(t)  is  a  single  valued  monotonically  increasing  function 
of  the  total  plastic  work  Wp(t)  such  that  an  increase  in  z(t) 
corresponds  to  an  increase  in  the  resistance  of  the  material 
to  plastic  flow.  This  implies  that  the  hardness  of  the  material 


is  taken  to  depend  on  the  amount  of  plastic  work  which  has  been 
done  on  the  material  from  a  reference  state.  In  (10)  ,  q  is 
a  material  parameter  describing  the  relative  amount  of 
hardening  that  is  isotropic,  and  ra(t)  are  taken,  following  [3] 
to  be  the  "direction  cosines"  of  the  rate  of  deformation  vector 

ra(t)  =  Da(t)/|D(t) |  (11) 

with 

|  D  |  2  =  D  •  D 

The  evolution  equations  (10)  for  the  state  variables 
which  control  the  hardness  of  the  material  indicate  that  part 
of  the  rate  of  change  is  isotropic  and  the  remaining  part 
varies  according  to  the  sign  and  orientation  of  the  current 
rate  of  deformation  vector.  When  q  =  1,  the  completely  iso¬ 
tropic  work-hardening  situation  is  obtained  in  which  all  the 
Aaa  are  equal  and  plastic  deformation  is  incompressible.  The 
case  q  <  0  is  admissible,  e.g.  cyclic  softening. 

Integration  of  the  evolution  equations  (10)  yields 

(o)  t  # 

Z  (t)  =  Z  +  e  /  qz  (x)  dx  + 
o 

t 

f a  sign [Da (t) ]  / (1-q) z (x) ra(x) dx  (12) 

o 

6 

where  Z(t)  =  £  Z  (t)e  is  a  6-dimensional  vector  described  by 
o=l  001  (o) 

the  base  vectors  ea,  and  Z  represents  the  initial  hardness  in 
each  of  the  six  coordinate  directions.  Equation  (12)  expresses 
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the  dependence  of  the  hardness  of  the  material  on  the  deformation 
path.  The  second  term  in  (12)  represents  the  total  isotropic 
hardening  up  to  time  t,  whereas  the  total  anisotropic  hardening 
used  to  characterize  the  multi-axidl  Bauschinger  effect  is  given 
by  the  last  term.  As  a  consequence  of  (12)  ,  the  components 
of  Z  at  the  current  time  t  have  different  values  which  gives 
induced  anisotropic  plastic  flow. 

Representation  of  the  anisotropic  plastic  flow  can  be 
completely  determined  by  specifying  Aaa  in  (8)  and  z  (t)  in  (10). 
Motivated  by  equations  relating  dislocation  velocities  and 
stresses,  Bodner  and  Partom  [1]  proposed  a  relation  which  can 
be  generalized  here  for  isothermal  anisotropic  hardening  in 
the  form 

V.  -  DoexP‘-iH1<ZL''3J2)nl/J21/2  (13> 

where  DQ  is  the  limiting  value  in  shear  of  the  strain  rate  of 
the  material  for  high  stresses  and  n  is  a  specific  parameter. 

The  function  z(t)  in  (10)  is  chosen  to  have  the  form 

4(t)  =  m(z1-zo)exp[-mWp(t)/zo]Wp(t)  (14) 

so  that 

z(t)  =  z1-(z1-z0)exp[-mWp(t)/zo]  (15) 

The  rate  of  the  plastic  work  is  given  by 


(16) 


6 

%  =  l  TaDa  +  trace  (oij)(Dl+D2+D3)/3 
c  a=l  J 

and  it  should  be  noticed  that  the  trace  is  an  invariant  quantity. 
In  (14) ,  zQ,  z^  and  m  are  appropriate  parameters  characterizing 
the  material. 

It  is  seen  from  (15)  that  an  increase  in  the  plastic 
work  Wp  corresponds  to  an  increase  in  z  which  implies  an 
increase  of  the  resistance  of  the  material  to  plastic  flow. 
Equation  (7)  and  (13)  would  lead  directly  to  isotropic  plastic 
flow  as  in  [1]  by  choosing  Zaa=  z  for  all  a.  For  the  general 
case  of  anisotropic  plastic  flow,  Zaa(t)  in  (13)  is  obtained  by 
substituting  z(t)  of  (15)  into  (12).  The  anisotropic  plastic 
representation  (12)  ,  (13)  and  (15)  requires  the  use  of  five 
material  constants  q,  n,  m,  zQ,  z^  which  must  be  determined 
experimentally  (DQ  is  actually  a  scale  factor) . 

Returning  to  equation  (7) ,  it  is  necessary  to  know  the 
orthogonal  transformation  which  diagonalizes  the  material 

A 

tensor  A^  of  (5) .  In  the  diagonal  form  A^  and  Za^  are  given 

with  respect  to  the  basis  vectors  e  .  If  the  material  is 

(o) 

initially  isotropic  Zaa(o)  =  Zaa=  zQ  for  all  a,  and  Za^=  Aa^=  0 
for  a^6  so  that  Aa^  is  in  diagonal  form.  Consequently,  any 
set  of  independent  vectors  efl  can  be  chosen  as  the  basis 
vectors.  As  the  deformation  proceeds,  the  Zaa  terms  determined 
from  (12)  become  different  from  each  other  and  give  rise  to 
anisotropic  plastic  flow  according  to  (13) . 


Formulation  of  the  Problem 


In  this  paper,  we  consider  the  dynamic  response  of  a 
sled)  -»<x^«»,  0<x2<H  under  plane  strain  conditions.  The  slab 
is  made  of  an  elastic  viscoplastic  anisotropic  work-hardening 
material  whose  constitutive  equations  were  described  in  the 
previous  section. 

The  slab  is  initially  at  rest  at  time  t<0.  The  surface 
Xj  =  H  is  assumed  to  be  kept  either  traction-free 

°12  ~  °22  =  0 '  -00  <  <  o°,  x2  =  H,  t  >  0  (17) 

or  rigidly  clamped 

u^  =  u2=0,  -°°  <  x^  <  °°,  x2=H,t>0  (18) 

Two  different  situations  for  the  slab  will  be  considered  in  this 
paper. 

(1)  Dynamic  indentation  by  a  rigid  punch: 

In  this  case/ the  surface  x2  =  0  of  the  slab  is  pressed 
at  t  =  0  by  a  rigid  indenter.  As  the  indenter  penetrates  the 
slab/  there  is  a  region  on  the  surface  x2  =  0  which  at  first 
has  traction-free  boundary  conditions  but  later  are  displacement 
conditions.  The  boundary  of  the  slab  conforms  to  the  geometry 
of  the  rigid  body  at  the  contact  region  which,  in  general,  is 
initially  unknown  and  must  be  determined  from  the  solution. 
Consequently,  the  appropriate  boundary  conditions  for  a  friction¬ 
less  contact  are 

_  _ _  _ _ _ _ _ .. _ - _ _ _ _ _ _ _ _ 
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°21  =  ° 


u2  =  f(xx,t) 


xll  1  X(t)  •  x2  ~  0'  t  >  0  (19) 


°21  =  0 


°22  =  0 


xx|  >  X(t) ,  x2  =  0,  t  >  0  (20) 


where  f(xlft)  is  the  prescribed  x2  displacement  imposed  over 
the  time-dependent  region  of  contact.  The  function  X(t) 
describes  the  boundary  of  the  moving  region,  which  is  a  priori 
unknown,  with  X(0)  =  0. 

(2)  Dynamic  excitation  by  time -dependent  normal  tractions: 


Here  we  assume  that  extended  normal  tractions  are  applied 
on  the  surface  x2  =  0  so  that 


°21  -  0 


°22  = 


-oo  <  Xl  <  °°,  x2  =  0,  t  >  0  (21) 


with  g  (x^,t)  describing  the  form  of  the  spatial  and  temporal 


excitation. 


The  stress-strain  relations  in  the  plane  strain  conditions 
for  which  e13  =  e23  =  e33  *  0  are  obtained  from  (2)  and  (1)  in 


the  form. 


aij  "  Xlell+e22]6ij+2yeij"Xekk)6ij"2yeij) 


The  displacements  u1(x1,x2,t)  and  u2(xlfx2,t)  are  governed  by 
the  dynamic  equations  of  motion.  These  follow  by  substituting 
(22)  into  the  equilibrium  equations: 


in  the  absence  of  body  forces ,  with  p  being  the  density  of  the 
material.  Substituting  (22)  into  (23)  leads  to  the  equations 
of  motion: 


"  (X+2y)ul,ll  +  yu1#22  +  (X+y)u2^12  -  \tWgl 


-2u(e{P>i  +  e^2) 


pii2  =  (X+2y)u2f22  +  yu2>n  +  (X+y)ulfl2  -  \e™2 


■2y(e^|l  +  ^22),2) 


For  isotropic  hardening,  plastic  deformation  is  incompressible, 
i.e.,  o  and  (22),  (24)  reduce  to  those  given  in  [2].  In 

the  following  section,  we  present  a  numerical  procedure  which 
enables  us  to  solve  the  above  mentioned  dynamic  problems. 
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Numerical  Treatment 

The  numerical  treatment  of  the  two-dimensional  dynamic 
problem  for  the  viscoplastic  slab,  which  was  formulated  in 
the  previous  section,  is  based  on  a  finite  difference 
procedure.  This  is  performed  by  introducing  a  net  of  mesh 
sizes  Ax^,  Axj  in  the  x^  and  x2  directions  respectively 
together  with  a  time  increment  At.  The  present  numerical 
treatment  is  a  generalization  of  the  procedure  given  in  [2]  for 
isotropic  plastic  flow. 

The  procedure  can  be  divided  into  three  parts.  In  the 
first  part,  the  displacement  components  at  interior  points 
within  the  slab  are  computed  from  the  equations  of  motion  (24) . 

In  the  second  part,  the  inelastic  variables  are  computed 
throughout  the  slab.  Finally,  the  third  part  consists 

of  the  computation  of  the  displacements  at  the  boundaries 
x2  =  0  and  x2  =  H. 

Part  1:  The  displacement  components  u1(x1,x2,t),  Ujtx^x^t)  at 
interior  points  within  the  slab  are  governed  by  the  dynamic 
equations  (24) .  These  are  solved  by  replacing  all  the 
derivatives  by  their  corresponding  central  difference  expressions 
obtaining  an  explicit  three-level  system  of  difference 
equations.  Thus,  it  is  possible  to  compute  the  displacements 
u^,  u2,  at  time  t  +  At  provided  their  values  at  the  previous 
and  current  steps  t  -  At  and  t,  respectively,  as  well  as  the 
values  of  plastic  strains  at  time  t,are  known  throughout 

the  medium.  The  form  of  this  difference  system  of  equations  is 
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essentially  similar  to  that  given  in  [4]  in  a  perfectly 
elastic  solid,  and  the  additional  inelastic  terms  add  no 
basic  difficulty.  This  system  is  of  a  second  order  accuracy, 
i.e.,  the  error  by  which  the  exact  solution  of  the  dynamic 
equations  (24)  does  not  satisfy  the  corresponding  difference 
equations  at  a  point  is  of  second  order  in  the  increments. 

Part  2 :  In  order  to  compute  the  inelastic  variables,  we 
denote  the  plastic  strains  by  E  (x.,x-,t)  such  that  E  =  D  with 

u  x  £  cx  m 

Dq  defined  in  (7) .  The  plastic  strain  and  the  plastic  work 
are  governed  by  (7)  and  (16)  respectively  with  A  given  in 

cx  cx 

(13) .  For  the  plane  strain  problem,  these  equations  provide 
a  system  of  five  equations  in  the  unknown  plastic  strains  and  the 
plastic  work  in  terms  of  the  variables  at  time  t  at  any  point 
in  the  region.  These  equations  can  be  integrated  to  yield 
the  Ea  and  at  time  t  +  At  at  a  point  of  the  slab  provided 
that  Aaa(x^,X2,t)  is  known  at  that  point. 

In  order  to  compute  Aaa(x^,x2»t)  at  a  P°int  at  time  t, 
the  hardness  Zaa(xi'x2't)  must  be  determined  [see  eqn.(13)]. 

This  is  achieved  by  employing  (10)  which  leads  approximately 
to 

Zaa(xl'x2't)  "  zaa(xl'x2't"At)  +  (x^x^t) 

+  (l-q)Az(x1,x2»t)r(x1,x2,t-At)  (25) 


where  Az  is  computed  from 


(x, ,x9,0)  =  z  assuming  that  the  material  is  initially 


Having  computed  Z 


to  determine  Z^tx^rXj/t)  .  This  is  then  employed  in  (13)  to 
obtain  ^aa ^xi fX2^  so  (7)  and  (16)  can  serve  to  provide 

the  values  cf  D  (x, ,x-,t)  and  W  (x, ,x~,t)  respectively. 


plastic  strains  and  plastic  work  at  time  t+At  by  using  the 
improved  Euler-Cauchy  method.  According  to  this  method 
E  and  W  at  time  t+At  are  tentatively  predicted  from 


E  (x. »x~ , t+At)  =  E  (x.,x0,t)  +  AtD  (x. ,x0 ,t) 


W  (x. ,x,,t+At)  =  W  (x.  f x. ft)  +  AtW  (x. » x_ , t) 


E  (x. tx~ , t+At)  =  E  (x. ,x- , t) 


+  At [D  (x, ,x~ ,t)  +  Dn(XwX,  t+At)l/2 
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Wp (xi /x2 / t+At)  =  W  (xlfx2,t) 

* 

+  tWp  (xi  ,x2  ft)  +  W  (x^x^t+At)  ]/2  (29) 

and  the  stars  over  D  and  Vf  indicate  that  these  are  evaluated 

a  p 

*  * 

using  Ea(x1,x2,t+At)  and  wp  ^  ,x2  ,t+At)  .  On  the  other  hand, 

*  *P 

the  evaluation  of  D  and  W  requires  the  determination  of 

*  p 

Aaa(xi»X2,t+At)  and  these  can  be  computed  after  obtaining 
+ 

Zact^xl'x2't+At^  *  Th^s  performed  in  the  same  manner  as  in  the 
computation  of  Z*a(x^,x2#t)  according  to  (25),  but  using  instead 
the  starred  quantities,  i.e., 

±  +  * 

Zaa*xl'x2't+At*  =  Zaa^xl'x2 +  qAz  (xi»x2't+At) 

+  (1-q)  Az  (x^  ,x2  ,  t+At)  r  (x^  ,x2  ,t)  (30) 

Consequently,  it  is  possible  to  compute  the  plastic  strains  and 
plastic  work  at  any  point  within  the  region  at  the  time  t+At 
when  all  the  field  variables  are  known  at  time  t-At  and  t. 

Part  3:  The  boundary  displacements  at  x2=0  are  determined 
according  to  the  relevant  boundary  conditions  [(19-20)  for  the 
indentation  problem  or  (21)  for  the  applied  tractions J .  For  the 
dynamic  indentation  of  the  slab  by  a  rigid  body,  the  boundary 
conditions  (19-20)  at  x2= 0  are  imposed  as  in  [2],  taking  into 
account  the  fact  that  the  region  of  contact  between  the  slab  and 
the  indenter  is,  in  general,  unknown  in  advance  but  must  be 
determined  as  part  of  the  solution.  This  was  achieved  by  an 
iterative  procedure  which  is  continued  until  all  the  requirements 
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of  the  indentation  problem  are  satisfied  simultaneously 
(see  [2]  for  more  details) . 

When  time-dependent  normal  tractions  are  applied  on 
the  surface  3^=0,  the  boundary  conditions  are  given  by  (21). 
These  are  imposed  by  approximating  the  derivatives  in  the 
x^ direction  by  central  difference  expressions  and  the 
derivatives  in  the  x2  direction  by  forward  difference 
expressions.  Consequently,  a  system  of  algebraic  equations 
in  the  unknown  displacements  at  the  boundary  x2=0 
obtained  at  every  time  step  as  follows : 

u1(iAx1,x2,t)  -  u1 (iAx1,x2+nAx2 ,t) 

-n{e [u2 ( (i+1) Ax1#x2 ,t)  -  u2 ( (i-1) Ax1#x2 ,t) ] 
-2Ax2E4 (x1,x2,t) }  =  0 

u2  (iAxlfx2  ,t)  -  u2  (iAXj^  ,x2+nAx2  ,t) 

-ri{e6  [ux  ( (i+1)  Axlfx2,t)  -  Uj^  ( (i— 1)  Axx  ,x2 ,t)  ] 

-Ax2  I261E2  (x^,x2,t)  +  6  (E1  (x1  ,x2,t)  +  EjtXj^Xjft) 
+E3  (x1,x2,t) )  +  gUj^t)]}  =  0 

where 

e  =  Ax2/(2Axx)  ,  6  -  A/(A+2y)  ,  6X  -  y/(A+2y)  , 
x2  *  o  ,  n  =  1. 


(31) 


Since  the  loading  under  consideration  is  symmetrical  about 


the  Xj  axis,  the  points  x^<0  need  not  be  considered, 

1 • e • ,  i=0 ,1,2 , .  . 

The  traction  free  conditions  (17)  at  the  surface  x2=H 
also  gives  the  system  of  equations  (31)  but  with  x2=H,  n=-l, 
g(x^,t)=0.  As  in  the  perfectly  elastic  problem  [5],  it  is 
convenient  and  efficient,  both  from  the  rate  of  convergence 
and  programming  viewpoints,  to  solve  the  above  system  of 
equations  by  the  Gauss-Seidel  iterative  procedure.  Details  of 
the  proof  of  the  convergence  of  the  process  are  given  in  [5] . 


I 


'  I 
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Applications 

The  proposed  method  has  been  applied  to  the  problems  of 
the  dynamic  indentation  of  a  viscoplastic  slab  by  a  long  rigid 
circular  cylinder,  and  to  that  of  a  slab  subjected  to  an 
extended  normal  load  rapidly  applied  on  the  slab  surface. 
Calculations  were  based  on  titanium  as  the  slab  material  for 
which  the  elastic  and  plastic  constants  are  given  in  [1] : 

X  =  0.9366  x  105  N/mm2,  u  =  0.44  x  105  N/mm2,  p  =  4.87  gm/cm3 , 
zQ  =  1150  N/mm2,  =  1400  N/mm2,  DQ  =  104  sec"1',  n  =  1  and 
m  =  100.  The  material  parameter  that  prescribes  the  relative 
amount  of  isotropic  hardening,  q,  was  found  in  (6]  to  be 
q  =  0.05  for  titanium.  This  was  obtained  under  uniaxial 
cyclic  stressing  conditions  using  the  current  direct  stress 
as  the  reference  for  the  relative  hardening  effects.  The 
value  q  =  0.05  was  used  in  these  calculations  although  it  is 
not  exactly  representative  since  the  hardening  "cosine  law" 
in  the  present  formulation  is  based  on  the  current  plastic 
deformation  rate  (11) .  Further  studies  on  this  subject  were 
made  in  relation  to  the  uniaxial  stress  problem  discussed 
later  in  the  paper. 

The  slab  thickness  was  taken  to  be  H  =  c  D  ”V 5  which 

o  o 

means  that  a  compressional  elastic  wave  whose  speed  in  the 

1/2 

material  is  cQ  -  [  (At2v0  /p]  will  propagate  the  distance  5H 
during  the  time  interval  D  Results  presented  in  this 

paper  were  obtained,  as  in  [2]* with  the  spatial  increments 
Ax^/H  =  AXj/H  =  0.1  and  the  time  step  cQAt/H  ■  0.05. 

(1)  The  dynamic  indentation  by  a  rigid  long  circular  cylinder 


21 


For  a  long  rigid  cylinder  indenting  the  slab,  the 
function  f(x^,t)  in  (19)  has  the  form 

f(x1#t)  =  p(t)  -  R  +  (R2-x12)1/2  (32) 

•  •  * 

In  this  equation,  p(t)  is  the  penetration  distance  of  the  indenter 
along  the  x2~axis,  i.e.,  p(t)  =  u2 (0 ,0, t) , and  R  is  the  radius 
of  the  cylinder.  This  is  assumed  to  be  much  larger  than  the 
contact  area  so  that  the  boundary  conditions  (19-20)  are 
applicable.  When  the  cylinder  is  pressed  into  the  slab  at  a 
constant  speed  V,  p(t)  =  Vt. 

In  Fig.  1,  plots  of  the  normal  and  tangential  stresses 
a22'  °11'  vertical  displacement  u2#  and  the  plastic  work  Wp 

are  given  versus  the  distance  x^  along  the  surface  x2=0  of  the 
slab  when  cQt/H=2.5  and,  5  with  =  0.05  and  R/H=5.  The  other 

surface  of  the  slab  is  kept  rigidly  clamped  so  that  (18)  holds. 

In  the  same  Figure,  the  plastic  work  along  the  plane  of 
symmetry  x^=0  within  the  slab  is  also  shown.  The  stresses  are 
normalized  with  respect  to  oQ  =  X+2y.  • 

These  results  are  compared  with  those  corresponding  to 
completely  isotropic  plastic  deformations,  i.e.  q«l.  This  com¬ 
parison  indicates  that  the  effect  of  the  anisotropic  hardening 
on  the  stress  field  and  the  shape  of  the  deformed  surface  of 
the  slab  is  not  significant  in  this  case. 

The  plots  of  the  plastic  work  Wp  give  an  indication  of 
the  amount  of  olastic  deformation  at  various  points  of  the  slab 
at  a  given  time.  The  anisotropic  hardening  case  is  seen  to  lead 
to  lower  values  of  Wp  compared  with  the  isotropic  one.  This 


does  not  seem  to  appreciably  influence  the  surface  stress  field 
nor  the  contact  area  or  the  force  required  to  press  the  cylinder 
into  the  slab  in  this  example. 

It  is  interesting  to  compare  the  amount  of  plastic 

3 

volume  change  E(x.,x-,t)  =  £  E  (x.,x_,t)  which  is  obtained  in 

1  a=l  a  1  1 

the  case  of  anisotropic  plastic  flow.  (When  the  flow  is 
isotropic  the  plastic  deformation  becomes  isochoric  at  all 
times  so  that  E(x^,x2,t)  =  0).  Table  I  presents  E(0,0,t)  at 
several  times  after  the  initiation  of  the  identation  process. 

The  Table  clearly  shows  a  monotonic  increase  with  indentation 
of  the  volume  change  due  to  plastic  die  formation  at  the  given 
point. 

The  present  indentation  problem  was  also  examined  for 
the  case  when  the  surface  x2=H  of  the  slab  is  kept  traction-free 
so  that  (17)  is  applicable.  Here  too,  the  results  for  the 
stress  field  show  similar  correspondence  between  the  aniso¬ 
tropic  and  isotropic  hardening.  Extensive  studies  of  the 
response  to  dynamic  indentation  of  a  slab  made  of  an  elastic- 
viscoplastic  isotropic  work-hardening  material  can  be  found 
in  [2]. 

(2)  Applied  time-dependent  loading: 

The  viscoplastic  slab  is  subjected  to  an  extended  normal 
time-dependent  traction  (21)  at  x2=0  of  the  form 


and  g  (x^,t)  =  0  for  |x^|>H.  In  Fig.  2,  the  resulting  stresses 
and  plastic  work  are  shown  at  x2=H/2  versus  x^  when  cQt/H  =  2.5 
and  5  with  go/cQ  -  0*1  and  c0tm/H  =  2.5  in  (33).  In  this  Figure 
we  present  also  the  plastic  work'  along  the  plane  of  symmetry 
x^=0  within  the  slab  at  the  same  times.  These  results  are 
compared  with  those  corresponding  to  the  isotropic  hardening 
case  (q=l) .  In  both  cases,  the  other  surface  of  the  slab 
x2=H  is  assumed  to  be  rigidly  clamped  (18) .  The  comparison 
shows  that  when  cQt/H  =  2.5,  the  resulting  stress  fields  are 
still  similar.  When  cQt/H  =  5,  however,  considerable  differences 
occur  in  the  stress  fields  between  the  two  cases  which  exhibits 
the  effect  of  anisotropic  hardening. 

In  Fig.  3 ,  we  present  the  results  for  the  same  two  cases 
as  Fig.  2  but  the  surface  x2=H  of  the  slab  is  kept  traction-free 
(17) .  Here  the  effect  of  the  anisotropic  hardening  is  less 
pronounced.  The  different  results  obtained  in  Figs.  2  and 
3  are  attributed  to  the  clamped  and  free  boundary  conditions 
which  affect  the  reflected  stress  waves  differently. 

Finally,  Table  II  gives  the  amount  of  plastic  volume 
change  at  the  point  x^=X2«0  at  several  times,  when  the  surface 
x2sH  of  the  slab  is  either  rigidly  clamped  or  traction-free. 


(3)  Uniaxial  deformation  due  to  dynamic  loading: 


In  order  to  obtain  better  understanding  of  the  effect 
of  anisotropic  hardening,  the  simplest  one-dimensional  problem 
of  uniaxial  dynamic  loading  of  a  thin  rod  of  length  H  in 
the  x2  direction  was  considered.  The  rod  was  assumed  to 
be  made  of  elastic-viscoplastic  anisotropic  work  - 
hardening  material  (titanium) .  In  this  case  the  only  non-zero 
stress  component  is  the  axial  stress  a22  so  that  (2)  reduces 
to 

^22  =  ^22^E 


where  E  =  y  (3X+2y)/ (X+y)  is  the  Young's  modulus  of  the  material. 
The  constitutive  equation  (22)  then  reduces  to 

a22  =  Ee22  "  Ee22 

with  e22  =  3u2/3x2,  so  that  (24)  take  the  form 


pu2  -  E(u2,22  -  e22^,2) 


(35) 


It  is  assumed  that  x2=0  is  loaded  in  the  form 


a22 


gL  sin(2nt/tm) 


< 


t  <  t_ 
—  m 


(36) 


L°  t  >  fcm 

with  g^  being  an  amplitude  factor,  and  the  other  end  is  kept 
stress-free,  i.e. 


This  one-dimensional  dynamic  problem  forms  a  special  case  of 


-1  - 3  - r - - 3l'“o  - - o  nr 

in  (36) ,  we  find  that  the  axial  stresses  when  q  =  0.05  are  the 


the  effect  of  anisotropic  hardening 


is  negligible  for  all  times.  It  turns  out  that  the  resulting 


stresses  saturate  in  this  case  so  that  the  effect  of  anisotropic 
hardening  would  therefore  be  relatively  small,  e.g.  17],  in 


accordance  with  our  observation.  On  the  other  hand,  choosing 


g^/aQ  =  0.005  and  cQtm/H  =  0.5  in  (36),  for  which  the 
amplitude  of  the  applied  stress  is  in  the  vicinity  of  the  yield 
stress,  the  effect  of  anisotropic  hardening  becomes  fairly 
pronounced,  Fig. 4. 


All  the  results  presented  in  this  paper  were  based  on  the 
assumption  of  [3]  that  the  anisotropic  change  in  resistance  to 
plastic  flow  is  in  the  direction  of  the  deformation  rate,  so 


that  the  "cosine  law"  (11)  for  the  rate  of  plastic  deformation 
is  applicable.  Another  possibility  is  to  use  a  "cosine  law" 
for  the  current  stress  which  means  that  r  (t)  in  (11)  would  be 


with  obvious  relations  between  a  and  i 


Note 


In  a  uniaxial  problem. 


that  (3J-) 


(38)  gives  r2  =  1  for  o22  >  0  and  r2  =  -1  for  o22  <0.  As 
previously  mentioned/  the  relative  amount  of  isotropic  hardening, 
q  in  (10) ,  was  determined  for  titanium  in  [6]  using  this  rule. 

It  would  be  interesting,  therefore,  to  compare  the  effects 
of  anisotropic  hardening  using  either  (38)  or  (11).  This 
comparison  also  appears  in  Fig.  4  which  exhibits  the  axial  stress 
and  plastic  work  versus  time  at  the  observation  point  H/2  at  the 
middle  of  the  rod  when  it  is  loaded  according  to  (36)  with 
g1/aQ  =  0.005  and  cotm/H  =0.5.  It  is  seen  that  these  two 
cosine  laws  do  not  lead  to  appreciable  differences  in  this  case. 

A  consequence  of  this  observation  is  that  the  value  of  q  chosen 
to  match  a  given  uniaxial  stress  response  should  not  be  sensitive 
to  the  choice  of  either  anisotropic  hardening  law.  Fig.  4  shows, 
however,  that  the  response  is  considerably  different  when  hardening 
is  assumed  to  be  isotropic  (q=l) . 

Finally,  Table  III  presents  the  amount  of  plastic  volume 
change  at  the  mid-point  of  the  rod  at  several  times  when  the 
"cosine  laws"  (11)  and  (38)  are  used.  The  Table  shows  that  the 
amount  of  volume  change  when  employing  the  stress  "cosine  law"  (38) 
is  about  three  times  larger  than  that  obtained  using  (11)  with 
the  same  value  of  q.  The  plastic  volume  change  is  therefore  a 
more  sensitive  indicator  of  the  difference  in  hardening  laws  than 
the  stress  response  for  the  uniaxial  stress  case.  Which  of  the 
two  references  for  anisotropic  hardening,  i.e.  stress  or  plastic 
strain  rate,  is  closer  to  the  physical  situation  is  not  obvious 
and  will  be  the  subject  of  further  investigations.  It  is  noted 
that  the  stress  rule  (38)  leads  to  simpler  numerical  work. 
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Conclusions 

Constitutive  equations  for  elastic-viscoplastic 
anisotropic  work -hardening  materials  proposed  by  Stouffer  and 
Bodner  [3]  have  been  implemented  to  study  the  effect  of  aniso¬ 
tropic  hardening  and  the  resulting  induced  anisotropic  plastic 
flow  in  several  two-dimensional  dynamic  problems.  A  basic 
assumption  in  the  formulation  [3]  is  that  anisotropic  hardening 
is  controlled  by  the  sign  and  orientation  of  the  rate  of 
plastic  deformation.  An  alternative  formulation  could  be 
adopted  by  taking  a  "cosine  law"  based  on  the  current  direct 
stress.  These  two  law  are  investigated  in  a  uniaxial  dynamic 
problem.  Further  studies  on  the  applicability  and  implications 
of  these  hardening  laws  are  under  investigation. 


TABLE  II 


time 

E (0#  Ort)  x  102 

c0t/H 

rigidly  clamped 
rear  surface 

traction-free 
rear  surface 

in 

• 

0 

0 

0 

1.0 

-0.21 

-0.21 

1.5 

-0.57 

-0.56 

2.0 

-0.91 

00 

in 

. 

0 

i 

2.5 

-1.14 

-0.44 

3.0 

-1.27 

1 

0 

• 

as 

as 

3.5 

-1.33 

l 

0 

• 

00 

0 

4.0 

-1.40 

-0.89 

4.5 

-1.48 

rH 

O 

• 

rH 

1 

5.0 

-1.59 

-1.19 

Variation  with  time  of  the  amount  of  volume  change 
due  to  plastic  deformation  at  the  point 
for  a  slab  loaded  at  x 2=0  when  the  other  surface  is 
rigidly  clamped  (18)  and  traction-free  (17) . 


: 


TABLE  III 


time 

c0t/H 

amount  of  volume  change 
due  to  plastic  deformation 

"strain  rate" 

"stress" 

cosine  law  (11) 

cosine  law  (38) 

0.5 

0 

0 

1.0 

-0.13  x  10"5 

-0.34  x  10"5 

1.5 

-0.26  x  10"3 

-0.69  x  10"3 

2.0 

-0.40  x  10-3 

-0.11  x  10-2 

2.5 

-0.81  x  10’3 

-0.22  x  10"2 

3.0 

-0.54  x  10"3 

-0.18  x  10"2 

3.5 

-0.53  x  10"3 

-0.18  x  10-2 

4.0 

-0.44  x  10-3 

-0.16  x  10"2 

4.5 

-0.35  x  10-3 

-0.14  x  10-2 

5.0 

-0.35  x  10”3 

-0.14  x  10“2 

Variation  with  time  of  the  amount  of  volume  change 
due  to  plastic  deformation  at  the  mid-point  of  the 
rod  based  on  the  "strain  rate"  cosine  law  (11)  and 
the  "stress"  cosine  law  (38)  for  anisotropic  hardening 
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FIGURE  CAPTIONS 


Figure  1: 


Figure  2: 


Figure  3: 


Figure  4 : 


Uniform  frictionless  indentation  of  an  elastic- 
viscoplastic  anisotropic  work -hardening  slab 

( - )  and  isotropic  work -hardening  slab 

( - ),  by  a  long  rigid  circular  cylinder  with  a 

speed  V/cQ  =  0.05.  The  plots  show  the  normal 
and  tangential  stresses ,  vertical  displacements 
and  plastic  work  versus  the  distance  along  the 
surface  x2=0,  and  also  the  plastic  work  along 
the  plane  of  symmetry  x1=0  within  the  slab,  at 
times  cQ  t/H  =  2.5  (top)  and  cQ  t/H  =  5  (bottom). 

Time -dependent  normal  loading  of  an  elastic- 
viscoplastic  anisotropic  work-hardening  slab 

( - )  and  isotropic  work -hardening  slab  ( - ) 

whose  rear  surfaces  are  rigidly  clamped.  The 
plots  show  the  stresses  and  plastic  work  versus 
the  distance  along  x2  =  H/2 ,  and  also  the  plastic 
work  along  the  plane  of  symmetry  x^=0  within  the 
slab,  at  times  cc  t/H  =  2.5  (top)  and  cQ  t/H  =  5 
(bottom) .  The  loading  function  is  given  by  (33) . 

Same  as  Figure  2  but  for  a  slab  whose  rear 
surface  is  traction-free. 

Axial  stress  and  plastic  work  versus  time  at  the 
mid-point  of  a  rod  made  of  elastic-viscoplastic 
anisotropic  work -hardening  material  whose  aniso¬ 
tropic  hardening  is  governed  by  the  "cosine  laws" 

on  (11)  ( . ),  and  on  a  (38)  ( - ),  and  for 

isotropic  work-hardening  material  ( - )  .  The 

loading  function  is  given  by  (36). 
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